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INTRODUCTION 

Since the first observations of neutron stars thirty years ago, they have affected 
many branches of physics. These extremely compact stars serve as natural physical 
laboratories for probing the properties of matter under extreme physical conditions. In 
particular, more than half of them possess magnetic fields B > 10 12 G. 

Despite their name, neutron stars consist not only of neutrons. They have a 
crust containing ionized iron, heavier elements, and exotic neutron- rich nuclei, 1 above 
which lie liquid and gaseous outer envelopes, which are thought to be composed of iron 
or lighter elements. 2 The atmosphere, that affects the spectrum of outgoing thermal 
radiation, likely consists of hydrogen, the most abundant element in the Universe, which 
might be brought to the star surface by fall-out of circumstellar medium. Neutral atoms 
can provide an appreciable contribution to the atmospheric opacity. 

Apart from the physics of neutron stars, quantum-mechanical calculations of strong- 
ly magnetized hydrogen atoms find application also in the physics of white dwarf stars 3 ' 4 
and in the solid state physics. 5 Because of this practical demand, hydrogen in strong 
magnetic fields has been well studied in the past two decades. 6 The peculiarity of the 
problem for neutron stars is that an atom cannot be considered abstractedly from its 
thermal motion. Indeed, neutron star atmospheres are hot (T ~ 10 5 — 10 6 K), so 
that typical kinetic energies of the atoms are non-negligible in comparison with typical 
binding energies. Taking the thermal motion into account is highly non-trivial, because 
an atom moving across magnetic field is equivalent to an atom placed in orthogonal 
electric and magnetic fields, so that the cylindrical symmetry is broken. 

At 7 3> 1, where 7 = %uj c /2 Ryd = 5/2.35 x 10 9 G»l and u c is the electron 
cyclotron frequency, the collective motion effects 7 ' 8 become especially pronounced. In 
particular, so-called decentered states (with the electron localized mostly in the "mag- 
netic well" aside from the Coulomb center) are likely to be populated even at the rela- 
tively high densities p > 1CT 2 g cm~ 3 typical of neutron star atmospheres. These exotic 
states have been predicted two decades ago by Burkova et al. 9 and studied recently by 
other authors. 10-12 
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Collective-motion effects on the usual "centered" states have been first consid- 
ered in frames of the theory of perturbation. 8,13 Non-perturbative results covering both 
centered and decentered states were subsequently presented for binding energies and 
wavefunctions, 14 ' 15 oscillator strengths, 15 spectral line shapes, 16 and photoionization 
cross sections. 17 None of these data, however, has been published in an easy-to-use 
form of tables or analytical expressions. 

In this contribution I propose approximate analytical expressions for the binding 
energies of the hydrogen atom arbitrarily moving in a magnetic field typical of neutron 
stars, 300 < 7 < 10 4 . This range is physically distinguished, since at weaker fields the 
spectrum is strongly complicated by multiple narrow anticrossings, 14 while the upper 
bound, 7 ~ 10 4 , corresponds to the onset of non- negligible relativistic effects. 18 



THEORETICAL FRAMEWORK 

Motion of the hydrogen atom in a magnetic field can be conveniently described by 
the pseudomomentum K = m p t p +m e r e — (e/c)Bx (r e — r p ), where the subscript i = e or 
i = p indicates electron or proton, respectively, r; = — (ik/m^Vi — (gj/mjc) A(rj) is the 
velocity operator, m 8 the mass, q p = —q e = e the charge, and A(r) the vector potential 
of the field. Gorkov and Dzyaloshinskii 19 have shown that in the representation in 
which all components of K have definite values, the relative motion can be described 
in terms of a one-particle Hamiltonian which depends on K. 

It is convenient to describe the centered states of the atom using the relative 
coordinate r^ ^ = r e — r p as independent variable and the axial gauge of the vector 
potential, A(r) = |B x r. For the decentered states, the "shifted" representation 19 
is more convenient. In the latter representation, the independent variable is rW = 
r e — r p — r c and the gauge is A(r) = |B x (r — [{m p — m e )/m i i] r c ). Here, r c = ^§2-B x K 
is the relative guiding center, and mn = m p + m e- 

Let us assume that B is directed along the z-axis. The z-component of the pseu- 
domomentum corresponding to the motion along the field yields the familiar term 
K^/2mn in the energy, while the transverse components Ki produce non-trivial ef- 
fects. Therefore we assume K z = and Ki = K hereafter. 

If there were no Coulomb attraction, then the transverse part of the wavefunction 
could be described by a Landau function $ ns (r| L L ' ) ), where is the projection of 
in the (:q/)-plane. The energy of the transverse excitation is 

E ns = [ n + {m e /m p )(n + s)}hu c , (1) 

where the zero-point and spin terms are disregarded. 

A wavefunction ip K of an atomic state \k) can be expanded over the complete set 
of the Landau functions 

^ ) (r (,) )=E*n.(rf ) )^ l ,(z) > (2) 

ns 

where 77 = or 1 indicates the conventional or shifted representation, respectively (a 
generalization to arbitrary 77 proved to be less useful 15 ). The one-dimensional functions 
g$). K are to be found numerically. The adiabatic approximation used in early works 9,19 
corresponds to retaining only one term in this expansion. 

A bound state can be numbered 15 as \k) = \n K , s K , is, K), where n K and s K relate 
to the leading term of the expansion (|j), and v enumerates longitudinal energy levels 

El !SK!V (K) = E K - E^ kSk (3) 
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Table 1. Parameters of the approximation (f|) at 10 1 < 7 < 10 4 . 



s 





1 


2 


3 


4 


5 


6 


7 


Pi 


15.55 


0.5332 


0.1707 


0.07924 


0.04696 


0.03075 


0.02142 


0.01589 


P2 


0.3780 


2.100 


4.150 


6.110 


7.640 


8.642 


9.286 


9.376 


PS 


2.727 


3.277 


3.838 


4.906 


5.787 


6.669 


7.421 


8.087 


Pi 


0.3034 


0.3092 


0.2945 


0.2748 


0.2579 


0.2431 


0.2312 


0.2209 


P5 


0.4380 


0.3784 


0.3472 


0.3157 


0.2977 


0.2843 


0.2750 


0.2682 



and controls the ^-parity: g$l. K (— z) = (—l) v g™. K (z). For the non-moving atom at 
7 > 1, the states v = are tightly bound in the Coulomb well, while the states v > 1 
are hydrogen- like, with binding energies below 1 Ryd. The states with n^0 belong to 
continuum at 7 > 0.2 and will not be considered here. 

At small pseudomomenta K, the states v = remain tightly bound and cen- 
tered, the mean electron-proton separation x being considerably smaller than r c (for 
the hydrogen-like states v > 1, however, x is close to r c at any K). The larger K, 
the greater is the distortion of the wavefunction towards r c , caused by the motion- 
induced electric field in the co-moving reference frame, until near some K c transition 
to the decentered state occurs, and the character of the motion totally changes. With 
further increasing K, the transverse velocity decreases and tends to zero, whereas the 
electron-proton separation increases and tends to r c . Thus, for the decentered states, 
the pseudomomentum characterizes electron-proton separation rather than velocity. 

At very large K the longitudinal functions become oscillator-like, corresponding 
to a wide, shallow parabolic potential well. 9 For a fixed u, this limit is reached at 
K ^> (u+^) 2 h/ ob, where as is the Bohr radius. Still at arbitrarily large K, there remain 
infinite number of bound states with high values of v whose longitudinal wavefunctions 
are governed by the Coulomb tail of the effective one- dimensional potential. 15 

The decentered states of the atom at K > K c ~ 10 2 au have relatively low binding 
energies and large quantum-mechanical sizes, I ~ K/l au; therefore they are expected 
to be destroyed by collisions with surrounding particles in the laboratory and in the 
white-dwarf atmospheres. In neutron-star atmospheres at 7 ~ 10 3 , however, the de- 
centered states may be significantly populated. This necessitates inclusion of the entire 
range of K below and above K c in the consideration. 



ANALYTICAL APPROXIMATIONS 

Binding Energies of the Non-Moving Hydrogen Atom 

Extensive tables of binding energies of the hydrogen atom at rest with respect 
to the magnetic field have been presented by Rosner et al. 20 and supplemented by 
other authors. 21-23 Recently, the accuracy ~ 10~ 12 Ryd has been achieved. 24 In the 
astrophysics, a lower accuracy is usually sufficient, and simple analytical estimates are 
often desirable. 

For this reason, we have constructed a fit to E^°\ where E^J V = — ^^(O), in a 
possibly widest range of 7. For the tightly-bound states, we have 

£& 7 ) = In (exp [(1 + s)~ 2 } + Pl [ln(l + P2V1)?) + Pz Ml + P^ 5 )? % d - ( 4 ) 

The parameters pi — P5 depend on s; they are listed in table |I[ This fit is accurate to 
within 0.1-1% at 7 = 10 _1 — 10 4 , and it also provides the correct limits at 7 — > 0. 
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Table 2. Parameters of <g) at 1 < 7 < 10 4 . 



V 


1 


2 


3 


4 


5 


6 


a u 
K 


0.785 
1.724 


0.578 
0.765 


0.901 
1.847 


0.631 
0.717 


0.970 
1.866 


0.660 
0.693 



For the hydrogen-like states, we use the asymptotic result 25 

(0) _ lRyd f N = {v + l)/2, S ~ T l for odd u, 

^sv ( N + S y> Wlieie I N = v/2, 5 ~ (I117)- 1 for even ja 1 J 

We have obtained the following fits to the quantum defect 5: for odd u, 5 = (a u +b uv fy+ 
0.0777) -1 , where a u « 1 and 6,, « 2; and for even v, 5 = [a v + 1.28 ln(l + ft^ 1 / 3 )] -1 , 
where a„ ~ | and 6^ ~ |. More accurate values of a u and &„ are given in table ||. At 
1 < 7 < 10 4 , errors of these approximations lie within ~ 10~ 3 . 

Binding Energies of the Moving Hydrogen Atom 

For the moving hydrogen atom in a strong magnetic field, the first analytical fit 
to E(K) has been published by Lai and Salpeter. 26 It is rather accurate for the ground 
state at K < K c but cannot be applied to excited or decentered states. 

We describe the longitudinal energy @ by the formula 

i, (K) (K) 

1 nsu\ )\ 1 + (K/K c y/<* 1 + (K c /Ky/ a ' { ' 

The two-term structure of (j^j is dictated by the necessity to describe the two physically 
distinct regions of K below and above K c . The parameter a has the meaning of the 
width of the transition region near K c in logarithmic scale of pseudomomenta. 

For the tightly-bound states, we parameterize the dependencies E^ (K) as follows: 

E$HK) = E&- * E&m= *™ , (7) 

Zm eS + qiK /ti 0s0 y'r 2 + rj + q 2 n 

where r* = r c /a-Q = K/(^ au), q\ and q 2 are dimensionless fitting parameters, and 
m e ff is the effective mass which is close to (but not necessarily coincident with) the 
transverse effective mass M^ sv obtained by the perturbation technique. At 7 > 300, we 
put qi = log 10 (7/300) if s = Oandgi = 0.5 otherwise, q 2 = 0.158 [ln((l+0.1s)7/215)] 2 / 5 , 
and a = 0.053 ln(7/150). For the effective mass, we have m e ff = mn [1 + (7/70)°°] , 
where c = 0.937+0.038S 1 - 58 and 70 = 6150(1 +0.0389s 3 / 2 )/[l+7.87s 3 / 2 ]. For the critical 



pseudomomentum, we have K c = [c\ + ln(l +7/7i)]y2mH^ ^ The parameters c\ and 
71 take on the values d = 0.81, 1.09, 1.18, 1.24 and 71 = (8.0, 3.25, 2.22, 1.25) x 10 4 for 
s = 0, 1, 2, 3, respectively. For s > 4, we put C\ = 0.93 + 0.08s and 71 = 6500. 

In figure |I] the above fitting formulae are compared with our numerical results 15 
and with the previous approximations. 26 The figure demonstrates that the present 
approximations are valid at any K from to infinity. Appreciable discrepancies occur 
only in narrow ranges of K near anticrossings. 

Now let us turn to the hydrogen-like states. Their binding energies are approxi- 
mated by the same formula (^]) but with slightly different expressions for E^ and E&>. 
For these states, M^ sv exceeds mn by orders of magnitude, and the perturbation method 
fails already at small K, 13 rendering the notion of the effective mass practically useless 
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Figure 1. Energy spectrum of the hydrogen atom moving across strong magnetic fields. Numerical 
values (dots) are compared with the present analytical approximations (full lines) and with 
previously published 26 ones (dashed lines). 



for the fitting. Thus we consider m eS as effectively infinite and put E^J U (K) = 
{y > f). Furthermore, the transition region is not well defined, and therefore K c and a 
lose their clear meaning and become mere fitting parameters. For odd states, we have, 

approximately, K c = (i/ 5/4 7 /170) a9 ^2m n E(°) and a = 0.66 + u/20. For even states, 



K c = v v /(7/600)m HJ E (0) and a = 0.66. 

The function E^(K) that describes the longitudinal energy at large K is now 



E®{K) = { (2 Ryd)- 1 r 2 + {2v + l)rf 2 + q 2 r* + l/E l 



3/2 



1/2 ^/^o)^ 1 



with q 2 = v 2 - 1 for odd u, and q 2 = v 2 + 2 V I 2 log 10 (7/300) for even v (at j > 300). The 
first and second terms in the square brackets ensure the correct asymptotic behavior. 15 



CONCLUDING REMARKS 

The analytical approximations for binding energies presented in this contribution 
depend continuously on two arguments — magnetic field strength and transverse pseu- 
domomentum. They are accurate, typically, within a few parts in 100-1000. The 
accuracy can be improved by almost an order of magnitude by optimizing the param- 
eters m e fr, K c , a, qi, q 2 in equations ©-(§) separately at each discrete value of 7. 
Tables of such optimized parameters have been obtained and will be published else- 
where, together with analytical approximations of geometrical sizes of various quantum- 
mechanical states of the moving atom and oscillator strengths of radiative transitions 
among them. The atomic sizes play important role in distribution of atoms over quan- 
tum states in a plasma and in their contribution to the plasma absorption coefficients. 
For example, a size of an atom may be used to evaluate effects of "unbounding" of 
electrons caused by random charge distribution in the plasma. For non-magnetized 
hydrogen plasma, an approximate treatment of these effects was revised recently; 27 for 
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the strong magnetic fields analogous work is under way. Eventually, the analytical 
estimates of ^-dependencies of the binding energies, atomic sizes, and transition rates 
help to generalize previously developed models of fully ionized atmospheres of magnetic 
neutron stars 28 to the more realistic case of partially ionized atmospheres. 
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